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Abstract. Many recent observational studies have concluded that planetary systems commonly 
exist in multiple-star systems. At least ~20 %, and presumably a larger fraction of the known 
extrasolar planetary systems are associated with one or more stellar companions. These stel- 
lar companions normally exist at large distances from the planetary systems (typical projected 
binary separations are on the orders lO'^ - 10* AU) and are often faint (ranging from F to T 
I^ ' spectral types). Yet, secular cyclic angular momentum exchange with these distant stellar com- 
panions can significantly alter the orbital configuration of the planets around the primaries. 
One of the most interesting and fairly common outcomes seen in numerical simulations is the 
opening of a large mutual inclination angle between the planetary orbits, forced by differential 
nodal precessions caused by the binary companion. The growth of the mutual inclination angle 
between planetary orbits induces additional large-amplitude eccentricity oscillations of the inner 
planet due to the quadrupole gravitational perturbation by the outer planet. This eccentricity 
1-^ ' oscillation may eventually result in the orbital decay of the inner planet through tidal friction, 
'r^' as previously proposed as Kozai migration or Kozai cycles with tidal friction (KCTF). This 
Q , orbital decay mechanism induced by the binary perturbation and subsequent tidal dissipation 
?H niay stand as an alternative formation channel for close-in extrasolar planets. 
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^ ' 1. Introduction 



As of June 2008, at le ast ■-^20 % of the ~240 stars hosting planets ar e components 
of multiple-star systems (JDesidera fc Barbierill2007t iRaghavan et al.l 120061 ). Most of the 



binaries that harbor planets are very wide, with typical sky-projected separations ranging 
from 10^ to IC* AU. Majority of the stellar companions to planetary systems are FGK- 
(^ • type stars, but there are also later-type companions, including 8 M-dwarfs and one T- 
OO . type brown dwarf. These distant, and often faint, companions, however, may still play a 
^^ ■ significant role in alterin g planetary orbits aroun d the primaries on a very long timescale. 
As first pointed out by IZucker &: Mazehl (|2002l ). the two populations of planets, those 



around single stars and those around components of binary systems, seem to exhibit 
, . , different orbital properties. The effect of secular binary perturbations is particularly 
C^ ' discernible in the distribution of planets' orbital eccentricities that are possibly excited 
by the Kozai mechanism (§ 12. ip . For example, a large fraction of planets with highly 
eccentric orbits (e > 0.7) are associated with stellar companions. Also, the population 
of planets in binaries shows a higher median eccentricity (0.25) compared to that of 
th e planets around s ingle stars (0.20). |jj Another curious trend that was first noticed 
by IZucker fc Mazehl ( 20021 ) and has remained robust still to date is the population of 



massive hot Jupiters (nip sini > 2Mj and P < 100 days) found exclusively among binary 

f observation al biases seem to influe nce the observed eccentricity distribution only marginally; 
see for example IShen fc Turnerl l|2008f ). 
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systems (there are a few recent exceptions such as X0-3b or HD 17156 b, for which the 
presence of stellar companions has not yet been ruled out). 

Interestingly, secular eccentricity-pumping mechanisms such as the Kozai mechanism 
are closely related to the formation of close-in (possibly transiting) planets. In addi- 
tion to the traditional planetary migration theory incorporating planet-disk interactions, 
there is a growing attention to an alternative mechanism involving circularization of 
initially eccentric orbits through tidal dissipation in the central stars. There are two 
major ways a planet can gain a sufficiently large orbital eccentricity that leads to the 
tidal circularization and orbital decay: mutual, impu lsive gravitational sca ttering be- 



tween planets upon dissipation of the nascent gas disk (jNagasawa et al.ll2008l). and adia- 



batic secular eccentricity gr owth induced by the Kozai mechanism (|Wu fc Murravll2003 : 



iFabrvckv fc Tremaind 120071 ). The latter scenario, often called Kozai migration or Kozai 



cycles with tidal friction (KCTF), has been adopted for single-planet systems in binaries 
and produced models that are consistent with the observations, given certain assump- 
tions in the initial conditions. In this study, we investigate the dynamical outcomes of 
multiple-planet systems secularly perturbed by stellar companions and find that Kozai 
migration process is significantly enhanced with the presence of additional planetary 
companions. 



2. Secular Three-body Interactions 

We focus on the secular evolution of hierarchical four-body systems that consist of two 
planets around a component of a wide binary system. It would help to break up the prob- 
lem into two secular three-body problems in order to provide analytical understandings 
of secular four-body interaction: (1) the secular evolution of two planets around the cen- 
tral star (nearly coplanar, not necessarily hierarchical), and the secular evolution of each 
planet due to the quadrupolc perturbation from the binary (not necessarily coplanar, 
but highly hierarchical) . The former has been formulated through the Laplace-Lagrange 
Secular Theory, and the latter through the Kozai mechanism. 

2.1. Secular Planet-Binary Interaction 

Binaries that contain planets are typically separated by more than 100 AU. In such wide 
binaries it is likely that the initial orbital angular momentum of the planet (i.e., the spin- 
axis orientation of the primary star) is not correlated to the orbital angular momentum 
of the binary (Hale 1994) . Thus, distribution of the initial inclination of the planetary 
orbit with respect to the binary plane is expected to be isotropic. The Kozai mechanism 
takes its effect only when the initial relative inclination is greater than the critical Kozai 
angle /koz = 39.2°. In such cases, planet's eccentricity (e), inclination (/), the argument 
of pericenter (cj), and the longitude of ascend ing node (VL \ starts to precess secularly 
all on similar timescales (called Kozai cycles; lKozailll962l ). The orbital eccentricity of 



a planet oscillates, even if it is initially close to zero, with a nearly constant amplitude 
(see for example the red-dashed curve in Figure [I]), t hat is to the lowest order a function 



only of the initial inclination /q ( Holman et al.lll997 ): 



v/l-5/3cos2/o. (2.1) 

The timescale of the eccentricity oscill ation (and the evolution of other orbital elements) 



is determined as ( Kiseleva et al.lll998l ) 
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where the subscripts and bin denote the primary and the companion stars, respectively. 
During Kozai cycles, a planetary orbit is slowly and repeatedly torqued by the binary, 
keeping the orbital energy and thus the semi-major axis constant. A cumulative effect of 
the torque results in the secular apsidal precession of the planet, i.e., Wp precesses on the 
timescale tkoz- In a typical planetary system there are other sources of perturbations that 
also induce secular apsidal precession of the planet, such as the general-relativistic effect 
caused by the central star, or secular perturbations from other planetary companions. 
If any of these extra perturbations induces apsidal precessions on the timescale smaller 
than tkoz, the torque from the binary companion secularly averages out to zero, and the 
planetary orbit will remain nearly circular. 

2.2. Secular Mutual Interaction between Planets 

Analytical solutions for the long-term evolution of a pair of planets around a star due 
to their mutual gravitational interaction is for mulated through the second-order Laplace- 
Lagr ange secular theory (L-L theory hereafter: iBrouwer fc Clemencell961tlMurrav fc DermottI 



1999f) . The L-L theory also yields secular evolutionary timescales of the planets. The two 



eigenfrequencies 5+ and g^ derived from the L-L theory correspond to the apsidal preces- 
sion timescales, Ti^^± « 2iT/g±. |jj This apsidal evolution timescale derived from the L-L 
theory can be compared to the Kozai timescale tKoz (Eq. [2.2| ) to determine whether the 
Kozai cycles are suppressed or not. If, for instance, a double-planet system in a binary 
has initial configuration such that Tt^^2 < tKoz,2, then the Kozai cycles on the outer planet 
will be entirely suppressed, and the binary companion will play no significant role in the 
secular evolution of the outer planet. 

The L-L theory also yields one eigenfrequency / which scales with the coupling strength 
of the planetary ascending nodes [J. This also translates to the nodal precession timescale 
of a double-planet system as tq « 2n/f. 



3. Non-rigid Evolution of Planetary Orbits in Binaries 

Newly formed planets emerging from the circumstellar gas disk are expected to orbit 
in a common plane (along the stellar equator). However, in the quadrupolar gr avitational 



poten tial of the binary, this coplanarity can be broken as the planets evolve (Ta keda et al. 



20081 ). If the nodal couphng of the planets (scales with rn) is strong compared to the 



nodal perturbation caused by the binary companion (scales with 2rKoz), planetary orbits 
can evolve rigidly, maintaining near-coplanarity. On the other hand, if the nodal coupling 
between the planets is weak compared to the nodal perturbation introduced by the stellar 
companion, a large mutual inclination angle starts to grow between the planetary orbits, 
thereby inducing Kozai cycles to the orbit of the inner planet. 

In Figure [U an initially zero mutual inclination (/12) between the planetary orbits 
grows periodically to as large as ^ 60°. As soon as the mutual inclination between the 
planets reaches 40° at i ~ 2Myr, the inner planet's orbital eccentricity ei starts to 
oscillate chaotically, reaching as high as ~0.95. This growth of /12 is solely due to the 
differential nodal pr ecession of the planetary orbits introduced by the stellar companion 
( Takeda et al.l I2OO8I ) . It is almost independent of how inclined initially the planetary 



orbits are with respect to the binary plane, unlike the three-body Kozai mechanism. In 

f In certain limits the apsidal evolution is described as a nonlinear combination of multiple 
eigenmodes, and this simple relation between t he eigenfrequency and the apsidal precession 
frequency does not hold. See (|Takeda et al.ll2008l ) 

X If an additional quadrupolar potential such as stellar oblateness is introduced, there will be 
two eigenfrequencies in the I-Q space. See lMurrav fc DermottI (|1999| ). 
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Figure 1. Numerically integrated evolution of a non-rigid double- planet system in a binary. 
The system evolved from the following initial conditions: for the inner planet (solid black curve), 
mi — 0.02Mj, oi = 0.9 AU and /i = 45°; for the outer planet (dashed red curve), 7712 = 0.6Mj, 
a2 = 12.7 AU and I2 = 45°; and for the binary companion, rribin ~ 0.7A'Iq, abin = 803 AU 
and Ebin ~ 0.78. The binary orbit remains at /bin = 0°. The dotted blue curve represents the 
relative inclination angle I12 between the two planets (scale shown on the right axis) . At around 
t = 14.5 Myr, the inner planet's pericenter enters inside 0.1 AU of the central star (horizontal 
dashed line in the bottom panel) and remains there for ~ 2 x lO^yr. 

four-body systems, the orbital eccentricity of the inner planet can still grow to a very 
large value even if the planet's initial inclination is below the critical Kozai angle. 

In an opposite case in which the nodal coupling of the planetary orbits is sufficiently 
strong compared to the nodal precession of the outer planet induced by the stellar com- 
panion, planetary orbits secularly rotate together, remaining on a common plane. In such 
dynamically-rigid configuration, /12 and ei remain small throughout the evolution. 



4. Kozai Migration and Formation of Close-in Planets in Binaries 

In dynamically non-rigid systems inner planet's eccentricity secularly oscillates due to 
the perturbation from the outer planet. During these Kozai cycles, without any dissipative 
mechanism the semi-major axis ai of the inner planet remains constant throughout the 
evolution. This means that at each eccentricity maximum, the planet's pericenter distance 
9i = ai(l ~ ei) reduces greatly. It is known that in the regime where qi < 0.1 AU, strong 
tidal interaction with the central star leads to damping of the orbital eccentricity and 
the semi-major axis, and the planetary orbit migr ates toward the cent r al star until it 
eventually halts at a tight, near ly circular orbit ( Kiseleva et al.l Il998l : IWu fc Murray 
2003t iFabrvckv fc Tremaind 120071 ) . For example, in the example in Figure [T] the inner 
planet is expected to start decaying its orbit around t « 15 Myr if tidal dissipation is 
included in the calculation. 

This so-called Kozai migration scenario has been adopted to pr oduce some dynamica l 
models that are consistent with observations (e.g., HD 80606b; IWu &: Murray! 120031 ) . 
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given relatively fine-tuned initial conditions. The difficulty is to achieve such a small 
pericenter distance during the Kozai cycles that the tidal dissipation significantly drains 
the orbital energy from the planet. Thus these dynamical models require very large initial 
relative inclination angles, almost 90°. Such a large angle is not unrealistic, but not likely 
to be applicable to many systems. 




Figure 2. The minimum pericenter distance in AU reached by the inner planet within 500 Myr 
of numerical integrations with different sets of /bin and a2 . Other orbital parameters and masses 
are fixed as mo = rribin = I.OMq, mi = O.lMj, m2 — l.OMj, ai = lAU, Obin = 750 AU, and 
fibin ~ 0.5. The darker-shaded area corresponds to the smaller pericenter distance achieved by 
the inner planet. Planets that reach pericenter approach smaller than ~0.1 AU are subject to 
tidal dissipation and expected to start migrating toward the central star. 



However, as explained in the previous section, in a more realistic planetary system 
containing more than one planet, such a tight requirement on the initial inclination 
angle is no longer necessary. Figure [5] shows the distribution of minimum pericenter 
distances achieved by the inner planet during 500 Myr of numerical integration with 
different sets of initial 02 and /bin- The parameter space is divided into two evolutionary 
classes depending on the location of the outer planet; in the region where 02 S 15 AU, the 
planetary system is dynamically rigid, and the inner planet's eccentricity remains small 
(pericenter remains large). The only exceptions in this region are when the planetary 
orbits are inclined from the binary by more than ^ 60°, in which case close encounter 
or possible orbital crossing between the planets induced by the large-amplitude Kozai 
cycles of the outer planet impulsively excites ei to large values. If, on the other hand, 
0-2 ^ 15 AU initially, then the system is dynamically non-rigid, and even if the planetary 
orbits and the binary plane are initially nearly coplanar, additional Kozai cycles will 
significantly excite the orbital eccentricity of the inner planet. 

Such a numerical coverage of 02 - /bin space shows similar morphology even if different 
initial binary parameters and planetary masses are selected; there is always a critical 
orbital separation a2,crit that separates dynamically-rigid and non-rigid evolutions. As 
illustrated in Figure [5] the formation of close-in planets through tidal circularization of 
eccentric orbits should be particularly efficient if the outer planet is initially located 
beyond a2,cnt- This boundary a2.crit can be analyt ically determined b y simple formulae, 
given the initial orbital parameters of the system ( Takeda et al.lbOOST ). 
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5. Properties of Close-in Planets Formed via Enhanced Kozai 
Migration 

As demonstrated in the previous chapter, a steUar companion can effectively stir up 
planetary systems from a large distance and excite the eccentricity of the innermost planet 
by propagating its perturbation inward through planets on a reasonable timescale. Once 
the planet's eccentricity becomes sufficiently large and the pericenter distance sufficiently 
small, the planet begins to migrate toward the central star until it settles on a tight, 
circular orbit. It is possible to identify close-in planets in binaries that formed via such 
an "induced Kozai migration" scenario from their final orbital properties. 

The leading hypothesis for the planet migration mechanism involves the loss of angular 
mom entum through interaction with the protoplanetary gas disk (e.g., iGoldreich fc Sari 
20031) . If the planet-disk interaction is the common migration channel for extrasolar plan- 



ets, majority of the close-in planets are expected to be found on orbits along the stellar 
equator. On the other hand, if a planet migrates via the combined effect of the Kozai 
cycles and tidal dissipation, the angular momentum of the final planetary orbit should 
be largely misaligned with respect to the stellar spin axis since the planet's inclination 
evolves during the Kozai cycles lQ. The sky-projected angle A of spin-orbit misalignment 
is now measurable for t ransiting planets through the Rossiter-McLaughlin (KM) effect 



e.g., IWinn et al.l 120051 ). Among the first 8 extrasolar planets with available RM mea- 



surements, two of them show large spin-orbit mis alignments; HP 17156 b (A — 62° ± 25° 



iNaritaet al1l2008h and X0-3b (A = 70° ± 15°; iHebrard et aP l2008h . Although these 



results still need to be verified by follow-up observations, once the spin-orbit misalign- 
ment is confirmed they will be the first strong candidate transiting systems that formed 
through Kozai migration. 

Another curious observational trend is the presence of massive close-in planets found 
predominantly among binaries. A possible indication of this trend is that massive plan- 
ets somehow do not migrate efficiently through planet-disk interaction. If that is true, 
then Kozai migration triggered by binary perturbations would stand as a strong alter- 
native to form close-in planets with a wide range of masses. Another possibly common 
migration channel, planet-planet scattering d uring the early stage of planet formation 
(IChatterjee et al.l 120081 : iNagasawa et al.l l2008l ) , also faces a significant challenge in mi- 
grating massive planets into tight orbits because it is usually the lighter planets that 
get scattered inward. The Kozai mechanism combined with tidal damping so far is the 
only process that does not depend strongly on planetary masses, consistent with the 
exclusive presence of massive ones in binaries. This is an other reason to susp ect that the 
two transiting pla nets HP 17156b (rrir, 



3.09+ ^I^M.T 



Gillonet al.ll2008f l and X0-3b 



(mp = ll-13Afj; Ijohns-Krull et al.ll2008l ; IWinn et al.ll2008l ) are the outcomes of Kozai 
migration, as both planets are fairly massive. 

The orbital migration process proposed in this study requires at least four or more 
bodies in a system, including two stars and two (or more) planets. Alternatively, for 
a handful of known close-in single planets in binaries, detections of secondary planets 
will significantly constrain the dynamical origin of those systems. Note that the induced 
Kozai migration is more common in hierarchical planetary systems where the second 
planet orbits beyond a critical radius a2,crit such that the planetary orbits evolve non- 
rigidly. Unfortunately, detection of additional planets beyond a few AU will take several 
years or decades of continuous radial- velocity monitoring. Once the presence of additional 



f Theoretically, orbital inclination can be also damped toward the stellar equator through 
tid al dis s ipatio n . However, it usually takes longer than the estimated age of the planetary system; 
seelHutI l|198lD : IWinn et all (|2005| ). 
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planetary companions in those systems are ruled in or out, however, the dynamical history 
of known close-in planets can be significantly constrained. 

This work was supported by NSF Grant AST-0507727. 
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